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ABSTRACT 

It has been recently claimed that there is no correlation between the distribution of 
galaxies and their peculiar velocities within the Local Volume (LV), namely a sphere 
of R = 7h _1 Mpc around the Local Group (LG). It has been then stated that this 
implies that either locally dark matter is not distributed in the same way as lumi- 
nous matter, or peculiar velocities are not due to fluctuations in mass. To test that 
statement a set of constrained N-body cosmological simulations, designed to repro- 
duce the main observed large scale structure, have been analyzed. The simulations 
were performed within the flat- A, open and flat matter only CDM cosmogonies. Two 
unconstrained simulations of the flat- A and open CDM models were performed for 
comparison. LG-like objects have been selected so as to mimic the real LG environ- 
ment. The local gravitational field due to all halos found within each LV is compared 
with the exact gravitational field induced by all matter in the simulation. We conclude 
that there is no correlation between the exact and the local gravitational field obtained 
by pairwise newtonian forces between halos. Moreover, the local gravitational field is 
uncorrelated with the peculiar velocities of halos. The exact gravitational field has a 
linear correlation with peculiar velocities but the proportionality constant relating the 
velocity with gravitational field falls below the prediction of the linear theory. Upon 
considering all matter inside the LVs, the exact and local gravitational accelerations 
show a much better correlation, but with a considerable scatter independent on the 
cosmological models. The main conclusion is that the lack of correlation between the 
local gravitation and the peculiar velocity fields around LG-like objects is naturally 
expected in the CDM cosmologies. 
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1 INTRODUCTION 

A key ingredient of the standard model of cosmology is 
that the large scale structure (LSS) has emerged out of an 
otherwise a homogenou s and isotropic universe via gravi- 
tational instability (e.g. IPeeblesI (Il980h ). One of the main 
consequences of gravitational instability is that the growth 
of structure induces a non- vanishing velocity field. The stan- 
dard model of cosmology relates the large scale mass den- 
sity and peculiar velocity fields. It further assumes that 
the observed galaxy distribution is closely related with the 
underlying matter density field, even if this relation is bi- 
ase d in some yet u nknown way. In an interesting recent pa- 
per IWhitinj (|2005h (W05) has challenged these basic ideas 
of the standard model by testing them against a study 
of the galaxy distribution and their peculiar velocities in 
the Local Volume (LV), defined as the sphere of radius 



R = 7h x Mpc (where h is the Hubble's constant in units of 
100km s _1 Mpc _1 ) centred on the Local Group (LG). W05's 
main conclusion is that 'Either dark matter is not distributed 
in the same way as luminous matter in this region, or pecu- 
liar velocities are not due to fluctuations in mass. ' This is 
a very important result and if it was true it would put the 
standard model of cosmology on very shaky foundations. 
The main goal of our paper is to examine W05 methodol- 
ogy and claims by analyzing the kinematics and dynamics of 
simulated LV sy terns in different cosmological models. As we 
will show in this paper, our numerical results do not support 
W05's claims. 

W05 carried out an analysis of the distribution and pe- 
culiar velocities of 149 galaxies in the LV. Using the high- 
quality data of these galaxies W05 has mapped the mass 
distribution within the LV, assuming it is traced by the 
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Figure 1. Projected dark matter distribution for the ACDM con- 
strained simulation. The image shows a projected slice of 10 h _1 
Mpc thick across the box centre in supergalactic coordinates X 
and Y. The Local Supercluster is the filamente crossing the box 
horizonally. The position of the Local Group and the Virgo and 
Fornax clusters are shown. 



galaxies, and calculated the gravitational field within the 
LV. The gravitational field has been calculated by summing 
over the pairwise Newtonian interaction for each galaxy and 
by weighting the galaxies by their K and B luminosity. W05's 
working assumption is that the peculiar velocity of galaxy 
should be aligned with the gravitational field it experiences 
and that the amplitude of the peculiar velocity and grav- 
ity fields should be linearly connected. The notion that 'the 
peculiar velocity field is linear with the gravitational field' 
is strictly valid in the linear regime of the gravitational in- 
stability in an expending universe. If such a linear relation 
had been confirmed by the high-quality data that the LV 
provides this would have validated the idea that the LSS is 
indeed induced by gravitational instability. No clear correla- 
tion was found between the velocity and local gravitational 
fields. This has led W05 to conclude that, at least in the 
LV, structure has not formed by means of the gravitational 
instability. 

The above conclusion should be tested carefully. The 
main goal of the present paper is to apply W05 analysis to 
a set of N-body constrained simulations (CSs) of the local 
universe. The CSs are designed to reproduce the gross fea- 
tures of the nearby universe, namely the main players of the 
nearby LSS, such as the Local Supercluster (LSC) and the 
Virgo , are imprinted onto the simulations (iKravtsov et al ' 
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initial conditions of the simulations as constrained realiza- 
tions of Gaussian fields, where actual observational data is 
used as constraints. The CSs provide an almost optimal lab- 
oratory for testing W05 algorithm, as they closely mimic the 
dynamics of the LV. CSs have been run in the framework 
of the benchmark flat-A CDM model as well as in the open 



CDM (OCDM) and the so-called flat matter only standard 
CDM (SCDM) models. In addition, W05 procedure has been 
tested against standard, non- constrained, simulations of the 
same models. 

This is a first in a series of papers that focuses on study- 
ing the nearby universe by means of N-body and hydrody- 
namical CSs. In particular issues concerning the coldness of 
the local flow, the mass distribution in the LV, the mass 
accretion history of the LG and the future of the nearby 
structure in a dark energy dominated universe are to be ad- 
dressed. The highlight of this project will be the full galaxy 
formation high resolution simulation of the LG. Our choice 
of the different models is dictated by our plan to study the 
nature of the velocity field in the LV in the CDM cosmogony 
and in particular its dependence on the dark energy and the 
dark matter. In this regard the SCDM model is considered 
to be an extreme case and is taken for reference. Forthcom- 
ing papers in this series will focus on the problem of the 
coldness of the local Hubble flow. 

The structure of the paper is as follows: a brief descrip- 
tion of the simulations is given in § [2j The selection of LG- 
like objects is described in §[3] and the analysis itself in § |4) 
The results of the analyses are given in § \5\ A final discussion 
of the results is given in § [6] 



2 SIMULATIONS 

W05 has analyzed one and only one particular patch of the 
universe, namely the LV that extends around the LG. The 
LG is not a unique or an unusual object in the universe, yet 
it has its own characteristics that affect the outcome of any 
dynamical test that would be applied to it. Our main goal 
is to check the validity of the W05 approach, in the context 
of the environment of the LG. This is to be achieved by 
applying the W05 analysis to LG-like objects found in ap- 
propriate N-body simulations. The key to a successful study 
is to properly select LG-like objects from the simulation and 
analyze the LV around these objects (see § [3]) . Yet the selec- 
tion of the LG-like objects is based on the properties of the 
LG itself and not on its environment within the LV. This has 
encouraged us to use CSs as a 'laboratory' for testing the 
W05 procedure in conditions very similar to the ones pre- 
vail in the LV and its immediate surrounding. In addition 
we have used unconstrained simulations for reference. The 
comparison of the results obtained in the constrained and 
unconstrained simulations will shed light on the question of 
whether the findings of W05 are naturally expected in CDM 
dominatted cosmogonies. 

The data used to constrain the initial conditions of 
the simulations is made of two kinds. The first data set is 
made of radial velocities of galaxies drawn from the MARK 
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HI dWillick et alJll997h , SBF 

iKarachentsev (|2005h catalogs 



Tonrv et all l200lh and the 
Peculiar velocities are less 
affected by non- linear effects an d are used as cons traints 
as if they were linear quantities (Zaroubi et al. 19 991). This 
follows the CSs pe rformed by Kravtso v et al.l |2002) and 
iKlvoin et all <l2Q03h . The other constraints are obtained 
from the catalog of nearby X -ray selected clusters of galaxies 
(|Reiprich fc BohrTn ger 2002). Given the virial parameters of 
a cluster and assuming the spherical top-hat model one can 
derive the linear overdensity of the cluster. The estimated 
linear overdensity is imposed on the mass scale of the cluster 
as a constraint. Different CSs with different random realiza- 
tions have been calculated and they all exhibit a clear and 
unambiguous LSC-like structure that dominates the entire 
simulation, much in the same way as in the actual universe 
in which the LSC dominates the nearby LSS. The simula- 
tions do vary with respect to the particular details of the 
LG-like object that is formed roughly in its actual position. 
All simulations used here are based on the same random 
number realization of the initial conditions. 

Five simulations have been performed so far. Three con- 
strained ones have been performed within the framework of 
the ACDM, OCDM and SCDM models. In addition, uncon- 
strained simulations of the ACDMand OCDM cosmologies 
have been performed for the sake of comparison and bench- 
marking. Table [T] presents the cosmological parameters of 
the different simulations. All simulations correspond to a 
periodic cubic box of 64 h _1 Mpc on a side. We made a 
random realization of the corresponding power spectrum for 
each cosmological model with a large number of particles 
(2048 3 ). At the same time, a constrained realization of the 
density field for the same power spectra was done in an eu- 
lerian mesh of 256 3 grid points. After fourier transforming 
both constrained and unconstrained density fields, we sub- 
stitute the fourier modes of the constrained field into the un- 
constrained one. Finally, we used Zeldovich approximation 
to compute the 3D displacement field for an initial redshift 
of z=60. Once we got the displacements for a refined mesh 
of 2048 3 grid points, we used it to estimate the initial condi- 
tions for a resampling of 256 3 dark matter particles in total. 
In this way, we are able to zoom into a particular area of the 
simulations and resimulate them with much higher resolu- 
tion, up to the maximum resolution possible (2048 3 ), having 
the same structures as in the low resolution simulations (see 
I Klvpin et all (|200lh for a more detailed information about 
zoomed simulation techniques). Thus, all the numerical ex- 
periments that are reported here have the same number of 
particles (256 3 ). This translates into a mass per particle of 
1.3 x 1O 7 /i _1 M for ACDMand OCDM and 4.3 x 1O 7 /i _1 M 
for SCDM simulations. 

We hav e used the pa rallel TREEPM N-body code 
GADGET2 (ISpringell <|2Q05h ) to run these simulations. An 
uniform mesh of 512 3 grid points was used to compute the 
long-range gravitational force by means of the Particle-Mesh 
algorithm. A constant comoving Plummer equivalent gravi- 
tational smoothing scale of 20 h~ x kpc was set at high red- 
shift and we changed it to 5 /i -1 kpc physical scale since 
z=3 till z=0. The number of timesteps to complete the evo- 
lution from z=60 till z=0 ranges from 5000 to 7000 depend- 
ing on the simulation. We employed a variety of parallel 
computer architectures (SGI-ALTIX, IBM-SP4, Opteron- 
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Table 2. Constrains used to find LG candidates following the 
Maccio et al. (2005) criterion. The circular velocity has been used 
for the mass cont rains. 



clusters) during the course of this work. Using 16 processors 
simultaneously, we completed one run in about 2 cpu days. 

In what follows, we will use the name ACDM, OCDM 
and SCDM for the simulations with Constrained Initial con- 
ditions in the different cosmological models. The names 
ACDMu and OCDMu will refer to the two different uncon- 
strained realizations in the ACDM and OCDM models re- 
spectively. As an example of how the simulations look like, 
we show in Figure [T] a projection of the dark matter distri- 
bution in the ACDM simulation box at z=0. 



3 SELECTION OF LG CANDIDATES 

Dark matter halos were found in simulations using two ob- 
ject finding methods: Th e Bou nd Density Maxima (BDM) 
algorithm (|Klypin et al.l (|l999h ) is based on finding local 
center of mass in spheres of variable radius starting from 
randomly sel ected part i cles in the simulation. The Amiga 
Halo Finder (jGill et al.l (|2004h ), on the contrary, finds local 
density maxima from an adaptive mesh hierachy. In both 
cases, an iterative procedure to find local centre of mass 
from density maxima is used. Particles that are not gravi- 
tatonally bound to the halo potential are also removed. We 
took halos composed by more than 100 dark matter par- 
ticles, which translates into a minimum mass per halo of 
Mm in = 1.3 x 10 9 /i _1 M o for the ACDM and OCDM simu- 
lations and Mmin = 4.3 x 1O 9 /i _1 M for the SCDM simu- 
lation. We indentified the same objects with both methods. 
For the work reported here we have used the halo catalogues 
obtained by the public available AMIGA halo finder code 
( [http://www.aip.de/People/AKnebe/AMIGA7] ). 

To identify LG candidates from the halo distribution, 
we selecte d those objects that ful fill t he strict requirement s 
as given in lGovernato et all (|l997h and lMaccio et al.l (|2QQ5h . 
They are summarized in Table [2] In brief, we searched for 
two halos similar to Milky Way and M31 galaxies, without 
neighbours with masses as high as any of the LG members 
and with a Virgo- like halo at an appropiate distance. A few 
tens of LG-like objects have been identified in each simula- 
tion (23 objects in ACDM, 34 in ACDMu, 41 in OCDM, 58 
in OCDMu and 37 in SCDM). One of the LG-like object for 
the ACDMsimulation which closely resembles the actual LG 
in terms of its mass and position, is presented in Figure [2] 
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Figure 2. Projected dark matter distribution around the best LG 
candidate in the ACDM simulation in supergalactic coordinates. 
The outer circle delimits the Local Volume and the inner circle 
represents the LG position 
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Figure 3. Plots of v vs. go and v vs. gL for the candidate of 
Figure E] 



4 ANALYSIS 

In order to study the dynamics of the LG-like objects found 
in the simulations, we have first computed the local Hubble 
flow in spheres of 7 h _1 Mpc around each candidate, and have 
also estimated the local overdensity within these spheres 
from the total mass inside them. All DM halos within the 
LV around each LG-like objects have been identified and the 
gravitational field acting on each halo has been calculated in 
two different ways. First, the local gravitational field field is 
calculated like in W05 by the summation over the pairwise 
newtonian interaction. Namely, the field acting on the i-th 
halo is given by 



-tot 



r0 2 + A 2 W 2 



(1) 



where, following W05, a softening parameter of A = 
1.2h _1 Mpc is introduced. Then we calculate also the 'true' 
gravitational field, namely the field calculated by the N-body 
code of the full mass distribution in the computational box. 



This is defined as the global gravitational field, (see the ap- 
pendix for further information). 

As in W05, the local gravitational field is decomposed 
into two terms, the contribution of a smooth background 
of matter and a fluctuating part given by the point mass 
distributions. Since we are interested in deviations from the 
average, we have to substract the linear term contributed by 
the background, g^, from the local gravitation field 



-tot 

g; 



(2) 



The calculation of the background term is done in two ways. 
First, following W05, the background solution is fitted by a 
linear term in r. Alternatively, the background solution is 
calculated by the exact solution for the unperturbed uni- 
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-G-^-Qrnpc^ 



(3) 



The two methods give virtually identical results and the 
fitting method has been used here, so as to be consistent with 
analysis of W05. We have also substracted the anisotropic 
background estimated from the tidal field (see § 15. 2[) . 

As we have mentioned earlier, the purpose of the present 
study is to compare gravitational accelerations and veloci- 
ties. To facilitate such a comparison the gravitational field 
is scaled by the linear theory prediction, 



gx 



2/(On,flA) _ 

3H Q 



(4) 



where x stands here for the local or global field. The velocity- 
gravity scaling factor is given by 



(5) 



(e.g. IPeeblesI (Il980h . lLahav et all (|l99lh ). Throughout the 
paper the gravitational fields will be represented by their 
scaled versions. 

Our analysis is exemplified by Figure [2] which shows 
one of the LG-like objects in the ACDM constrained sim- 
ulation. The figure shows the matter distribution in a box 
of 20h _1 Mpc centered on the simulated LG and projected 
on the Supergalactic plane. In Figure [3] two scatter plots 
show the scatter of the local and global gravitational accel- 
erations, vs. the peculiar velocity of all DM halos in this 
particular LV. In what follows, all quantities shown corre- 
spond to projections along the line of sight with the observer 
located in the center of mass of the LG candidates. 



5 RESULTS 

5.1 Relation between local and global 
accelerations 

Figure [4] shows scatter plots of local vs. global gravitational 
fields for some of the Local Volume candidates extracted 
from simulations. In order to study the relation between 
both fields, a linear fit to the gL vs. gG distribution has 
been made for all the LG candidates found in the different 
simulations. 

Figure [5] shows the least square fit slopes and correlation 
coefficients of the local to the global gravitational fields as a 
function of the overdensidity Sp/p measured within each of 
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Figure 4. Local vs. global scaled acceleration for one of the LG candidates in the constrained ACDM, OCDM and SCDM simulations 
respectively. 
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Figure 5. Slopes (first row) and Pearson's correlation coefficient (second row) of the linear fit of gi, vs. go as a function of overdensity 
calculated within the simulated Local Volumes for all candidates found in each simulation. The median and the average slope and 
correlation coefficient with their one sigma errors for all LGs in each simulation are also shown. 



the simulated LV's. We find that the local and global field 
are uncorrelated, showing an extremely small correlation co- 
efficient and a mean and median slope of roughly 0.1 rather 
than unity. 

We have tested the effect of taking only haloes with high 
mass (more than 10 11 h _1 M©). The results are very similar 
to the previous ones, because the lightest haloes do not con- 
tribute much to the gravitational field. The overdensity was 
calculated by computing all matter inside the local volumes. 
When using only mass in halos to estimate the overdensity, 
results do not change at all. 

In the preceding analysis the local gravitational field is 
calculated by summing over the pairs of DM halos, assum- 
ing they carry all the mass in the LV and that their mass 
is spread out over a relatively large area described by the 
A smoothing parameter. The poor correlation between the 
local and global gravitational fields has led us to relax the 
assumption that the mass is traced by the DM halos and 
calculate the local field contributed by all dark matter par- 



ticles in the LV. This is presented in Figure [6] which shows a 
very tight correlation between gz, and gc for one particular 
LG-like object in the ACDM CS. A better correlation with 
a slope very close to unity are obtained when we include all 
the particles that constitute the inter-halo medium, rather 
than by the intra-halo particles only, as can be seen in Figure 
[71 Moreover, the only gravitational smoothing done in this 
case correspond to that included in the simulation, which 
is of the order of kiloparsec scale, contrary to the strong 
smoothing of 1.2 Mpc used in W05 analysis and in our pre- 
vious estimates. The horizontal branch around the origin 
observed in Figure [6] is mainly due to particles not bound 
to halos. These particles are more affected by the external 
field than particles bounds to halos, and therefore for these 
the local acceleration is much smaller that the global one. 

Obviously, the correlation between the local and global 
acceleration improves as the LV increases. This has been 
checked for the LG-like object of Figure [6] Assuming a LV 
of a 30 h _1 Mpc radius the correlation improves to S = 0.47 
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Figure 6. Local vs. global gravitational accelerations computed 
taking into account all the particles within the LV for the same 
candidate in the ACDM simulation as in Figure [4] The small 
points represent the individual dark matter particles. Thick solid 
points correspond to halos in which accelerations have been com- 
puted by the average of all particles inside them. The straight 
line shows the equivalence between both accelerations. 
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5.2 Accelerations and velocities 

Figures [8] and [9] present the peculiar velocities as a func- 
tion of the scaled global and local gravitational accelera- 
tions around some selected LG-like objects in the different 
simulations. The slopes and the correlation coefficients of 
the linear fits for the relation between the peculiar veloci- 
ties and the gravitational accelerations acting on halos in the 
LV around all LG-like objects are given in Figures [TOl and [TT1 
The least square fit analysis confirms the visual impression 
of Figures [8] and [9] The peculiar velocity of halos in the LV 
of LG-like objects is clearly correlated with the global accel- 
eration but with a non-negligible scatter, with a mean and 
median correlation coefficient around 0.5 for the constrained 
simulations and somewhat lower for the unconstrained ones 
fFigure fTO]) . The mean and the median slope of the linear re- 
lation between the velocity and scaled global acceleration is 
in the range of 0.4 to 0.6. The distribution of the fitted slopes 
shows some dependence on the mean overdensity (Sp/p) in 
the LV. The width of the distribution decreases with the 
overdensity. At low Sp/p's the slope ranges from roughly 0.2 
to almost unity but at Sp/p > 1 the slope shows a narrow 
scatter around its mean value of « 0.5. Upon scaling of the 
gravitational field the linear theory predicts the slop to be 
unity. It follows that the amplitude of the peculiar velocities 
is smaller than what is expected by the linear theory. 

The peculiar velocities do not show any correlation with 
the local accelerations (Figures l9l and [TT]) . Their linear fit 
yields extremely low correlation coefficients and the fitted 
slope has no meaning. One should recall here that this lack 
of correlation has already been anticipated from the lack of 
correlation we found between the local and global accelera- 
tions (§EQ). 

The calculation of the local gravitational acceleration 
neglects the contribution of the tidal field, which is induced 
by the inhomogenous matter distribution outside of the LV. 
Obviously, for objects like the LG, with the Virgo cluster 
located just outside of the LV, the tidal field cannot be ne- 
glected. The local gravity- velocity correlation should be im- 
proved by adding the tidal field into the fitting procedure. 
Indeed, we follow W05 and extend the fitting procedure to 
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Figure 7. Slopes (first row) and Pearson's correlation coefficients 
(second row) for the vs. go least square fits for all the LG can- 
didates found in the ACDM and OCDM simulations as a func- 
tion of overdensity. All particles within simulated LVs have been 
taken in the computation of the local accelerations. A gravita- 
tional smoothing parameter of 5 kpc is assumed to compute 
the newtonian pairwise forces between particles. 



and r 2 = 0.58. However this is still considerably worse than 
the S — 0.89 and r 2 = 0.91 obtained by considering all the 
particles within the 7 h -1 Mpc LV (Figure [6]). 

The lack of correlation between the local and global 
gravitational accelerations, strongly leads to expect the ab- 
sence of correlation between the local gravitational field and 
the peculiar velocities. This is clearly confirmed by the anal- 
ysis presented in § 15.21 



where vo is an unknown vector and H is a symmetrical 
tensor with six unknown quatities. The nine free parameters 
are found by minimizing the scatter. Note that in the limit 
of the linear theory and for a small LV vo should be set to 
zero and the symmetric tensor H should be traceless. As 
the above assumptions do not hold for the LG we allow for 
a finite vo and for H to have a trace. 

Figure IT2l shows that indeed adding the vo and H terms 
improves somewhat the correlation between the local grav- 
ity and peculiar velocities. Yet, they are weakly correlated 
with r 2 « (0.1 - 0.2) for most LG-like objects. The minor 
improvement is not surprising. The Virgo cluster is located 
at a distance of about 10h _1 Mpc and the LV is denned by 
a sphere of radius of 7h -1 Mpc. Hence, modeling the tidal 
field by a spatial linear expansion, as is implicitly assumed 
in Eq. [6] constitutes a poor fit of the tidal field. 
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Figure 8. Peculiar velocity vs. global scaled acceleration of halos inside the LV for the same candidates as in Figure |4] 
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Figure 9. Peculiar velocity vs. local scaled acceleration for halos within the LV for the same candidates as in Figure 0] 
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6 DISCUSSION 

We have clearly demonstrated in this paper that the W05's 
analysis is not expected to yield a simple and clear corre- 
lation between the gravitational field calculated from the 
mass distribution within the local volume and the peculiar 
velocities of halos within that volume. This result invali- 
dates W05's basic assumption that the peculiar velocity field 
traces the gravitational field in a simple fashion, and there- 
fore the statements concerning the role of the gravitational 
instability in the LV are not valid. 

Setting aside issues concerning the practical limitations 
posed by observations of the LV and the uncertainty in es- 
timating the dynamical mass of luminous galaxies there are 
three main theoretical reasons why W05's analysis fails. The 
most obvious one is that the gravitational field is assumed 
to be traced by the galaxies, treated as point-like particles. 
Figs. [4] and [5] show the poor correlation between the actual 
gravitational field and the local field induced by the galax- 
ies. Now, this poor correlation is mostly due to the sampling 
of the field by the galaxies, as manifested by Fig. [6] which 
shows a clear correlation between the global field and the 
local field that is induced by all dark matter particles in the 
LV. 

The other reasons for the breakdown of the simple grav- 
ity - velocity relation of the linear theory are both related to 
the tidal field, hence the sheer of the velocity field. The sim- 
pler reason is that in solving the Poisson equation one should 



not neglect the homogenous solution, namely the tidal field. 
Now, in principle this can be easily corrected by adding a 
(spatial) linear term to the gravitational field that scales 
with the traceless shear tensor (Eq. [6] and Figure [12]). Thus 
by adding six free parameters to the fitting procedure one 
might be able to account for the tidal field. However, the 
size of the LV is such that the spatial linear expansion of 
the tidal field would fail and lead to an incorrect estimation 
of the gravitational field. For the LV centered on the LG, 
the tidal effect of the LSC cannot be represented by a linear 
term. 

The other reason for the inadequacy of the linear the- 
ory is more subtle. It has been shown that in the quasi- 
linear regime the growth of the dens ity contrast depends 
on the magnitude of the shear tensor (lHoffman|[l9 86, 1989; 



Zaroubi fe Hoffman! Il993l: I van de Wevgaert fc Babul 



1994 



Bertschinger &: Jainl 19941 ). The shear dependence intro- 
duces a non-local term in the equations that govern the 
growth of structure in the quasi-linear regime. Indeed, Fig. 
[8] shows a tight linear relation between the peculiar veloci- 
ties and the scaled global gravitation field. Yet, the constant 
of proportionality is less than unity for all the LG-like ob- 
jects in all the simulation (except of one single object) as 
is pr edicted by the linear theory prediction (in agreement 
with iHofTmanl (| 19891 )1. This behavior implies that under the 
optimal conditions of a full knowledge of the gravitational 
field a linear relation between the gravity and velocity field 
is expected within the LV around LG-like objects. However, 
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Figure 10. Same as Figure [5] but for the peculiar velocity vs. global scaled aceleration fits. 
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Figure 11. The same as Figure [5] but for the peculiar velocity vs. local scaled acceleration fit. 



using the slope of the relation as a way of measuring 
would underestimate its true value. 



W05 attempted to find a simple linear relation between 
the peculiar velocities and local gravitational field of galaxies 
in the LV. A careful analysis of the dynamics within the LV 
around LG-like objects identified in constrained simulations 
of the local universe and in unconstrained simulations in 
flat-ACDM, OCDM and flat-matter only CDM cosmologies 
shows that a lack of correlation is to be expected. Hence we 
cannot support the claim that 'either dark matter is not dis- 
tributed in the same way as luminous matter in this region, 
or peculiar velocities are not due to fluctuations in mass. ' 
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8 APPENDIX: THE GADGET ESTIMATION 
OF GRAVITATIONAL ACCELERATIONS 

Throughout the paper the global and local gravitational 
accelerations have been compared. The local accelerations 
within the LV were calculated by summing the Newtonian 
pairwise field. This can be repeated for the global field by 
summing over all the particles of the simulation. This di- 
rect sum is computationally very costly and does not take 
into account the contribution of the infinite periodic boxes. 
We used instead the particle accelerations calculated by the 
TREE-PM algorithm by the GADGET code. The relation 
between the GADGET acceleration with gG is as follows: 

The physical r and comoving x coordinates are related 
by: r = ax. The global gravitational field equals the physical 
acceleration of an object r = g. 

Now, the GADGET code provides an acceleration-like 
term defined as: 



ld_ 

a dt 



(a • v p ) = v p + H • v p 



(7) 



where v p is the peculiar velocity and a is the expansion scale 
factor. 

It follows that 



r = fp + r-. 

a 

Recalling that 



3 



one gets at the end 



gG 



(8) 



(9) 



This is the value we used for the global acceleration of each 
dark matter particle. The total acceleration of halos was 
computed by averaging this quantity over the entire number 
of particles belonging to each halo. 
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